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Abstract. The quantum correlations, including entanglement and discord with its 
geometric measure, and classical correlation are studied for a bipartite partition of 
a open or closed quantum system. It is found that the purity of the initial state 
plays an important role in the dynamics of quantum and classical correlations. In the 
dephasing model, the quantum correlations loss and the classical correlation gain are 
instantaneously happen. While, the purity of the initial state destroys the quantum 
correlations which is resulted by the unitary interaction. Therefore, with the purity 
parameter, a particular region in which there is no state have quantum correlations 
can be determined. 
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1. Introduction 

The development of the quantum information technology stimulates a deep study of 
the properties of quantum correlations inherent in a quantum system. In particular, 
there is a problem of identification of those quantum correlations which are responsible 
for the advantages of the quantum computations in comparison with the classical ones. 
Therefore, lots of interest have been devoted to the definition and understanding of 
correlations in quantum systems in the last two decades. Specially, the definition and 
study of quantum and classical correlations in quantum systems. It is well known 
that the total correlation in a bipartite quantum system can be measured by quantum 
mutual information PQ, which may be divided into classical and quantum parts [2-5]. 
The quantum part is called quantum discord(QD) which is originally introduced by pE]. 
Recently, it has been aware of the fact that quantum discord is a more general concept 
to measure quantum correlation than quantum entanglement (QE) since there is a 
nonzero quantum discord in some separable mixed states [I] . The dynamics of quantum 
discord and entanglement has been recently compared under the same conditions when 
entanglement dynamic undergoes a sudden death [6j [7J El El [TO]. It was shown that 
quantum discord presents an instantaneous disappearance at some time points in non- 
Markovian regime [6], and exponential decay in Markovian regime [7]. Interestingly, 
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it has been proven both theoretically and experimentally that such states provide 
computational speedup compared to classical states in some quantum computation 
models [HI [12]. In these contexts, quantum discord could be a new resource for quantum 
computation. 

The calculation of quantum discord is based on numerical maximization procedure, 
it does not guarantee exact results and there are few analytical expressions including 
special cases [T3| [T4"] . To avoid this difficulty, geometric measure of quantum discord 
(GMQD) is introduced by Ref.[15j. which measures the quantum correlations through 
the minimum Hilbert-Schmidt distance between the given state and zero discord state. 

Because of the unavoidable interaction between a quantum system and its 
environment, understanding the dynamics of quantum and classical correlations (CC) is 
an interesting line of research [14-17]. The influence of Markovian [7] and non-Markovian 
[T7] environment on the dynamics of QD and GMQD. They showed that both GMQD 
and QD die asymptotically with entanglement sudden death, and the discontinuity in 
the decay rate of GMQD does not always imply the discontinuity in the decay rate 
of QD. Also, they observed that even when QD vanishes at discrete times, GMQD 
disappears but not instantly. But in our work the dephasing environment, in which 
energy transfer from the system to the environment does not occur, is considered. 
Some work has been devoted to this issue [201 I2H 122], in which, the authors show 
that disentanglement is dependent on the initial condition and temperature of the 
environment. This type of the environment leads to phenomenon entanglement sudden 
death, i.e., the entanglement can decrease to zero abruptly and remains zero for a finite 
time [231 I2H 125]. Entanglement sudden death has been experimentally observed in an 
implementation using twin photons [26], and atomic ensembles |27j . 

In this paper, one considers a dephasing model, in which two qubits embed into a 
multi-mode quantized field and the interaction between the two qubits is also considered. 
Therefore, the quantum correlations via quantum discord and its geometric measure 
(GMQD) is compared with both quantum entanglement and classical correlation. 

2. Measures of correlations 

2.1. Quantum discord 

To quantify the quantum correlations of a bipartite system, no matter whether it is 
separable or entangled, one can use the quantum discord [21 H]. Quantum discord 
measures all nonclassical correlations and defined as the difference between total 
correlation and the classical correlation with the following expression 

D(p AB )=l(p AB )-Q(p AB ), (1) 

which quantifies the quantum correlations in p AB and can be further dis- 
tributed into entanglement and quantum dissonance (quantum correlations excluding 
entanglement) [28] . Here the total correlation between two subsystems A and B of a 
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bipartite quantum system p is measured by quantum mutual information, 

l(p AB )=S(p A )+S(p B )-S(p AB ), (2) 

where S(p AB ) = Tr(p AB log p AB ) is the Von Neumann entropy, p A = Ttb(p ab ) 
and p B = Tta(p ab ) are the reduced density operators of the subsystems A and B, 
respectively The measure of classical correlation is introduced implicitly in Ref.[l] 
and interpreted explicitly in the Ref.[2]. The classical correlation between the two 
subsystems A and B is given by 

Q(p AB ) = max[5(/) - 5>SfoO], (3) 
{n fc } k 

where {Hk} is a complete set of projectors to measure the subsystem B, and pk = 
Tr B [(I A (g) Ii k )pAB{I A ® n fc )]/p fc is the state of the subsystem A after the measurement 
resulting in outcome k with the probability pk = TrAB[(I A <8> ^-k)p AB (I A <S> !!&)], and 
I A denotes the identity operator for the subsystem A. Here, maximizing the quantity 
represents the most gained information about the system A as a result of the perfect 
measurement {n^}. It can be shown that quantum discord is zero for states with 
only classical correlations and nonzero for states with quantum correlations. Note 
that discord is not a symmetric quantity, i.e., its amount depends on the measurement 
performed on the subsytem A or B [T5] . 



2.2. Geometric measure of quantum discord 

The geometric measure of quantum discord quantifies the quantum correlation through 
the nearest Hilbert-Schmidt distance between the given state and the zero discord state 
[T5l fT6] , which is given by 

D A =^\\p AB - X \\\ (4) 

where 5* denotes the set of zero discord states and ||v4|| 2 = Tr(A^A) is the square of 
Hilbert-Schmidt norm of Hermitian operators. The subscript A of D 9 A implies that the 
measurement is taken on the system A. A state x on H A ® H B is of zero discord if and 
only if it is a classical-quantum state [29], which can be represented as 

2 

X = J2Pk\ k )( k \ ® Pk, 
k=l 

where {pk} is a probability distribution, \k) is an arbitrary orthonormal basis for H A 
and pk is a set of arbitrary states (density operators) on H B . An easily computable 
exact expression for the geometric measure of quantum discord is obtained by Ref. [15] 
for a two qubit system, which can be described as follows. Consider a two-qubit state 
p expressed in its Bloch representation as 

p AB = 1 {I A ^ jB + £ Xt{(Ti ®JB + JA^ m(Ti) 
4 i=l 

' X , ' : J°' °J.- ( 5 ) 

ij=l 
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where {a{\ are the usual Pauli spin matrices. The components of the local Bloch vector 
are x { = Tr(p AB (ai <g> I)) and y t = Tr(p AB (I ® tr<)). i2y = Tr(p AB (ai <g> aj)) are the 
components of the correlation matrix[15]. Therefore, its geometric measure of quantum 
discord is given by 

^ = ^(PII 2 + IMI 2 -^), (6) 

x = (x±, X2, %3) T , R is the matrix with elements and k max is the largest eigenvalue 
of the matrix K = xx 1 - + RR T . 

2.3. Entanglement via Negativity 

Here, one uses the negativity [3U| to measure the entanglement, i.e., the negative 
eigenvalues of the partial transposition of p AB are used to measure the entanglement of 
the qubits system. Therefore, the negativity of a state p AB is defined as 

JV(p) = max(0, -2 $>,•), (7) 

3 

where pj is the negative eigenvalue of (p AB (t)) TB , and Tb denotes the partial transpose 
with respect to the second system. 

3. Quantum and classical correlations in two-qubit models 

In this section one tries to present, by examining two examples, an physical 
interpretation for the relation between quantum entanglement, quantum and classical 
correlations. The first model consists of two non-interaction qubits couple with the 
same quantized field under the rotating-wave approximation. Another is the dephasing 
model, in which two qubits embed into a multi-mode quantized field and the interaction 
between the two qubits is also considered. 

3.1. Two non-interaction qubits couple with the field 

Here, one considers two qubits coupled with a single-mode cavity field, which is in 
the case of the exact resonance with the qubits. One of the pioneering potential 
applications on this Hamiltonian in the context of quantum information is " Cooper 
pair box" (qubits) [3"T]. i.e., the coupled system of two Cooper pair box (artificial atoms) 
and photons stored in the resonator (cavity mode). The cavity is sustaining non- 
decaying single mode field in its thermal state along with the mode structure of the 
electromagnetic field. The interaction picture Hamiltonian with the rotating wave 
approximation is given by 

# = AX)(a+|0) fc <l|+a|l)*(0|), (8) 
fc=i 

where a) and a are the creation and annihilation operators for the cavity mode, |0)fe 
and denote to the ground and excited states of the fc-th qubit, respectively. For the 
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whole system (system+field) , the evolution of the whole system is characterized by the 
interaction between the two-qubits system and single-mode cavity. However from the 
point of the qubits, the energy transfer between the qubits and the field happens, which 
is described by the relaxation term $\ 0)^(1 1 and the backaction term a | l)fc(0|. One 
assumes that the two qubits are initially prepared in Werner states, which is defined by 

p AB (P)=p\<p)(<p\ + \{l-p)I, (9) 

where \tp) = sin 9\ 11) + cos#|00) and p is a real number which indicates the purity 
of initial state, J is a 4 x 4 identity matrix. But the cavity field is initially prepared 
in the vacuum state, i.e., p F {0) = |0)(0|. Then the initial density operation for the 
whole qubits-field system is: p ABF (0) = p AB (0) (g> |0)(0|. By using the above initial 
states, the density matrix of the qubits-field system with the interaction (jSJ) evolves to 
p(t) = U(t)p ABF (0)U^(t), where the time evolution operator U(t) = exp(-iHt). The 
reduced density matrix, p AB (t), of two qubits is calculated by tracing out the cavity 
field variables. Therefore, p AB (t) is given by 

p AB (t) = ai|ll)(ll| + a 2 |00)(00| + a 3 (|ll}(00| 

+ |00)(ll|) + a 4 (|10)(10| + |01)(01|) 

+ a 5 (|10)(01| + |01}(10|) (10) 

with the abbreviation 

1 — p p 2 • 2 a 

ai — — h - (2 + cos wt) sm 9, 

4 9 

P 2 i\ 2p , \ 2 2 /i 

a2 = — h p cos 9 + -g- (1 — cos wtj sin 9, 

p 

as = - (2 + cos wt) sin 9 cos 9, 

1 — p p . 2 .2/1 1 — p 

a± = — 1 — sm wt sm 9, a§ = a± , 

4 6 4 

where w = 2.4495A. The eigenvalues of the density matrix p AB (t) are given by:Ai j2 = 



a 4 ± a 5 , and A3 t 4 = \[{a\ + a 2 ) ± — a2) 2 — 4a§]. After some straightforward 
calculation, the reduced density matrices associated with the above states is given by 

p A (t)= P B (t) 

= ( ai + a 4 )|ll)(ll| + (a 2 + a 4 )|00)(00|. (11) 

One noted that, the reduced density matrices of the qubits are represented in diagonal 
matrices. Therefore, these the states p A (t) and p B {t) are classical states. 

For the case of the two qubits coupled with a single-mode cavity field, the results 
are given in FigsH^i,b. In these figures, one reports the dynamics of D 9 A , D(p AB ), N(p) 
and Q(p) as a functions of the time At for different values of the purity of initial state 
(namely p = 0.5,1.0) with 9 — j. It is worth noting to mention that, because the 
matrix of the initial state of the two qubits is not a diagonal matrix, this state is not 
classical state. Therefore, its quantum correlation have non zero value at Xt = 0.0. 
For p = 1.0, these measures instantaneously oscillate and reach their maximum (at 
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Figure 1. Time evolutions of the quantum discord (dash plots), the geometric measure 
of QD (sold plots), the negativity (dotted plots) and the classical correlation (dash-dot 
plots) for p = 1 in (a) and p = 0.5 in (b) for 9 = \. 

Xt = ^f-, n — 0, 1, 2, ...) and minimum values (at At = ( 2n+1 ) 7r ) a t the same time points 
(see Fig{Th). Because the function cos cut is a periodical function on the scaled time 
with period — , D 9 A , D(p), N(p) and Q(p) evolve periodically with respect to the scaled 
time with period — (see FigsOJi,b). A rather counterintuitive feature of the QE is that 
it may exceed the measures of the GMQD and QD (see FigQJ,). So one can say that 
GMQD and QD are more general than QE. 

For p = 0.5 (see Fig{Tb), one can observe that the phenomenon of entanglement 
death occurs, but this phenomenon does not occur for GMQD and QD even when 
the purity p is small. Because the entanglement undergoes sudden death while the 
correlations are long lived, QE is not greater than QD, GMQD and CC for some 
time. Also one can see that one value of QE corresponds to many values of QD 
and GMQD, meaning that the states in possession of the same entanglement give 
different correlations. This means that, there are correlations(quantum and classical) 
in the intervals of entanglement death. Therefore, the entanglement is not the only 
part of quantum correlations. This agrees with Ref.[29], which showed that absence of 
entanglement does not imply classicality. From Fig. la, one can see that the intervals of 
vanishing negativity disappear when p = 1. Therefore, the entanglement sudden death 
is completely sensitive to the purity of the initial-state. 

3.2. Dephasing two interaction qubits by a multimode quantized field 

Here, one considers a dephasing model of two qubits embed into a multimode quantized 
field and the interaction between the two qubits is also considered. The dephasing 
channel case is important situation for the open systems, in which there is no energy 
transfer between the system and environment. The Hamiltonian of this case can be 



Quantum discord and its geometric measure with death entanglement 



7 




Figure 2. Time evolutions of the quantum discord (dash plots), the geometric measure 
of QD (sold plots), the negativity (dotted plots) and the classical correlation (dash-dot 
plots) for 7 = 0.0A in (a) and 7 = 1.0A in (b)for 8 = and p = 0.5. 



written as [33] 

+ Z)w fc S£5 fc + lk(0-%(TB + °B°A)( h i + h), (12) 
k 

where u is the qubit transition frequency and ujk the frequency of the fc-th field, 
the coupling constant between the two qubits is A but 7& is system-reservoir coupling 
constant (dephasing channel parameter). 

By using the same previous initial states but with \ip) = sin 6 1 ! 10) + cos 6*101), 
the density matrix of the dephasing qubits-field system is given by: U(t)p ABF (Q)U\t). 
Therefore, the reduced density matrix of the two qubits is given by 

p AB {t) = ai(|ll)(ll| + |00)(00|) + a 2 |10)(10| 

+ a 3 |01)<01| + « 4 |10)(01| + a*|01)(10|, (13) 

with the abbreviation 

a 1 = {l-p)/4, a 2 = (l+p)/4-/3, 
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Figure 3. The same as in FigHbut for 6) = (2, |j) in (a) and (J, 6) = (2, §) in 
(b) 



a 3 = (1 +p)/4 + /3, a 4 = -sin 26* - 

/3 = -L<j cos 2$ cos 2\t, 

L d = exp[-4^(— ) 2 (1 - cosa; fc )]- 

k U k 

Where Ld is the decoherence factor which leads to damping of the off-diagonal terms. 
One can get the eigenvalues of the density matrix p(t) as 

1 



Ai l2 = ai, X 3)4 = -[l+p±yJ(l+p) 2 -16\a A \ 2 ] (14) 

The reduced density matrices associated with the above states are given by 

p A (t) = ( ai + a 2 )|ll)(ll| + («! + a 3 )|00)(00|, (15) 
p B (t) = (ai + a 3 )|ll)(ll| + (ai + a 2 )|00)(00|. (16) 

The results of the case of the dephasing model are given in Figsj2]-??. In Figl21 GMQD, 
QD, QE and CC are plotted as a functions of the time At for different values of system- 
reservoir coupling parameter (namely 7/A = 0.0, 1,2) with 9 = ^ and p = 0.5. From 
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Fig|2K one can easily find the common features of the dynamics of GMQD, QD, QE and 
CC. The previous measures of quantum and classical correlations present instantaneous 
oscillations and reach their extreme values at the same time points. Because the 
entanglement undergoes sudden death for 7^0, while D A , D(p) and Q(p) are long 
lived. N(p) is not greater than D A , D(p) and Q(p) for small intervals and QE attains 
constant values during an intervals while the quantum correlations (GMQD and QD) 
vary in these intervals. 

In Figsj2bl3h,b, one examines the effect of system-reservoir coupling 7 on the 
dynamics of the previous measures of quantum and classical correlations with 6 = 
It is clear that D A , D(p) and N(p) decrease with increasing of 7, and they have zero 
values. Precisely, the dephasing parameter 7 leads to exponentially decay for maximum 
values of the D A , D(p) and N(p) to zero value, while Q(p) exponentially evolves to 
its asymptotic value. It is interesting to note that the larger the value of 7 is, the 
more rapidly D 9 A , D(p) and N(p) reach its asymptotic values of zero. This means that, 
quantum correlations, including entanglement and discord with its geometric measure 
die asymptotically with large values of 7. On the contrary, the classical correlation 
increases with increasing 7, and it have nonzero values for 7 > 0. The figures show that 
the classical correlation approaches an almost steady state for large values of 7 faster 
than that for small values of 7. Finally, after a very long time, the classical correlation 
loses it oscillations and asymptotically reaches its steady state, i.e., the final state of 
the qubits reaches a classical state. One can say that, the quantum states with large 7 
are rapidly transformed into classical states, i.e., the processes of quantum correlation 
loss and classical correlation gain are instantaneously happen. 

To investigate the influence of 6 — | (mixed state) with a large value of 7 = 2A on 
the correlations see FigfJUb. From this figure, one notes that both QD and CC reach 
their asymptotic values and approximately have the same behavior while GMQD and QE 
have the same behavior. This shows that the mixedness of the initial states affects on all 
the previous measures in a similar way and it inhibits them from going into zero. In 
Figs.??,??, one examines the effect of the purity of the initial states on the dynamics 
of the previous measures with 6 = | and 7 = 0.8A. It is clear that they decrease 
with decreasing of the purity parameter p. When the purity p is zero, all the measures 
vanish, i.e, the mixedness of the initial states have the same effect on all measures. One 
sees that the influence of purity leads to: the amplitudes of the local maxima of D A , 
-D(p), N(p) and Q(p) have exponential decay with decreasing the parameter p. When 
quantum correlations measures quite vanishes, the states p AB finally go into a classical 
state and its quantum correlation is lost completely. This means that, after a particular 
time, the purity destroys the quantum correlations of the qubits which is resulted by the 
unitary interaction. Therefore, in the presence of purity one can determines a particular 
region, in which, there is no state have quantum correlations. 
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4. Conclusions 

The dynamics of quantum correlations, including entanglement and discord with its 

geometric measure, and classical correlation in two-qubit models are introduced for a 

open or closed quantum system. It is found that the dynamics of GMQD, QD, QE 

and CC differ. Where, quantum discord and its geometric measure are exist in the 

region where the entanglement is zero, which is a strong signature for the presence of 

non classical correlations. System-reservoir coupling leads to: GMQD, QD and QE 

die asymptotically with larger system-reservoir coupling parameter. Also, processes 

of quantum correlation loss and classical correlation gain are instantaneously happen. 

It is found that the purity of the initial states destroys the quantum correlations by 

exponential decay. Therefore, in presence of the purity, one can determines a particular 

region in which there is no state have quantum correlations. 
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